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ABSTRACT 

Simulations of the formation of galaxies, as well as ionisation models used to interpret 
observations of quasar absorption lines, generally either assume ionisation equilibrium 
or ignore the presence of the extra-galactic background (EGB) radiation. We intro- 
duce a method to compute the non-equilibrium ionisation and cooling of diffuse gas 
exposed to the EGB. Our method iterates the ionisation states of the 11 elements 
that dominate the cooling (H, He, C, N, O, Ne, Si, Mg, S, Ca, & Fe) and uses tabu- 
lated ion-by-ion cooling and photo-heating efficiencies to update the temperature of 
the gas. Our reaction network includes radiative and di-electric recombination, col- 
lisional ionisation, photo-ionisation, Auger ionisation, and charge transfer. We verify 
that our method reproduces published results for collisional equilibrium, collisional 
non-equilibrium, and photo-ionised equilibrium. Non-equilibrium effects can become 
very important in cooling gas, particularly below 10 6 K. Photo-ionisation and non- 
equilibrium effects both tend to boost the degree of ionisation and to reduce cooling 
efficiencies. The effect of the EGB is larger for lower densities (i.e. higher ionisation 
parameters). Hence, photo-ionisation affects (equilibrium and non-equilibrium) cool- 
ing more under isochoric than under isobaric conditions. Non-equilibrium effects are 
smaller in the presence of the EGB and are thus overestimated when using collisional- 
only processes. The inclusion of the EGB alters the observational diagnostics of diffuse, 
metal-enriched gas (e.g. metal absorption lines probed in quasar sight lines) even more 
significantly than the cooling efficiencies. We argue that the cooling efficiency should 
be considered if ionisation models are used to infer physical conditions from observed 
line ratios, as the a priori probability of observing gas is lower if its cooling time is 
shorter. We provide on-line tables of ionisation fractions and cooling efficiencies, as 
well as other data, for equilibrium and non-equilibrium scenarios, and both with and 
without an EGB. Cooling efficiencies and diagnostics of the physical state of diffuse gas 
can become highly inaccurate if ionisation equilibrium is assumed or if the existence 
of the ionising background is ignored. 

Key words: atomic processes; plasmas; galaxies: formation; intergalactic medium; 
quasars: absorption lines; cosmology: theory; 



1 INTRODUCTION 

The rates at which diffuse gases cool radiatively are central 
to numerous baryonic processes in astrophysics. In the con- 
text of galaxy formation, how gas cools in the intergalactic 
medium (IGM) and galactic haloes critically determines if, 
how, and when galaxies receive their gas, which provides the 
main source of fuel for star formation. 

The ejection of the nucleosynthetic products of star 
formation (metal species) into diffuse regions complicates 
the picture, as these sub-dominant heavier atomic species 
are very efficient coolants, owing to their more numerous 



electrons and line transitions. The inclusion of radiative 
metal cooling in gas enriched to solar metal abundances 
( Z©) or even 0.1 Z© significantly reduces cooling times 



(e.e. lOooen 


leimer fe Davd 120061; IChoi fe Nasa 


mine] 2009: 


Schave et al. 


2010: Ivan de Voort et al. 2012: 


Haas et al.l 


201%. 



Diffuse regions outside of galaxies are photo-ionised 
and photo-heated by the extra-galactic background (EGB) 
originating from a combination of UV and X-ray emitting 
sources including ac tive galactic nuclei (AGN) a nd star- 
forming galaxies fe.g. lHaardt fe Madaulfl 9961 . 1 20121 ). Photo- 
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ionisation of hydrogen, helium, and metal species reduces 
cooling efficiencies, because fewer bound electrons provide 
fewer line transitions for radiative cooling (|Efstathioulll992l : 
IWiersma et~aT1l2009aL hereafter W09). W09 showed net cool- 
ing efficiencies declined by an order of magnitude at typical 
densities and temperatures of the shock- heated IGM. In- 
cluding the effect of photo-ionisation on metal cooling ef- 
ficiencies can therefore significantly affect the appearance 
of m etal absorption in quasar absorption line (QAL) spec - 
tra (|Tepper-Garcia et al.1 l201ll : lOppenheimer et al.1 l2012h 
and alter the dynamics of the accreting ga s onto galaxies 
(jSmith et al.ll201ll : lOppenheimer et al.ll2012h . 

Most hydrodynamic simulations and semi-analytic 
models assume that gas is in ionisation equilibrium, 
but this assumption can break down if the ionisation 
or recombination time of a species is long compared to 
other timescales, including the dynamical, Hubble, or 
cooling timescales. For example, iGnat &; Sternberg! ([2007, 
hereafter GS07) calculated that when the timescale of 
recombination exceeded that for cooling, a "recombination 
lag" occurs where t he gas is over-ionised relative to the 
equilibrium case (e.g. iKafatos 1973; Sha piro fe Moorelll976l; 
Schmutzler &; Tscharnuterl 1993k ISutherland &; Dopita 



data used in this work at http: / / noneq.strw.leidenuniv.nl) 
Supplementary figures are also available at this website. 



19931 : ICen fe FandfeOOd : lYoshikawa fe Sasakill2QQ6l : IVasiliev 
2013). By following the ionisation state of gas without 
any external radiation (i.e. collisional processes only), they 
find fewer bound electrons for metal ion species with long 
recombination times and that this results in less efficient 
net cooling and cooling functions of temperature with less 
sharp peaks. 

In this paper, we develop a new method to follow 
non-equilibrium ionisation and cooling in the presence of 
a photo- ionising EGB, focusing on temperatures ^ 10 4 K. 
By self-consistently tracking the ionisation state of primor- 
dial and metal species, we can accurately calculate the non- 
equilibrium cooling rates of metal-enriched gas. In §2, we in- 
troduce and test the method t hat we have int egrated into the 
Gadget-3 (last described in ISpringell l2005h hydrodynami- 
cal simulations code for hydrogen, helium, electrons, and the 
important metal coolants (carbon, nitrogen, oxygen, neon, 
magnesium, silicon, sulphur, calcium, and iron). To demon- 
strate the effects and importance of non-equilibrium pro- 
cesses, we use a stand-alone, "single-particle" version of the 
code to follow the ionisation and temperature evolution of a 
parcel of gas independent of any hydrodynamics, but still in 
the presence of an EGB. We diagnose how non-equilibrium 
effects diverge from the equilibrium case for enriched cool- 
ing gases without photo-ionisation in §3 and with photo- 
ionisation in §4. Relevant astrophysical applications include 
metal-enriched materials ejected into the IGM and the cir- 
cumgalactic medium (CGM) by superwinds, which simula- 
tions s how regularly re- accrete onto galaxies at late time s 
[z < 1: lOppenheimer et alJlioiol ; Ivan de Voort et al.ll2012h . 
We show that the combination of non-equilibrium effects 
and photo-ionisation can retard accretion and significantly 
alter the observational diagnostics of cooling gas that may 
be observed via metal-line absorption with instruments such 
as the Cosmic Origins Spectrograph (COS) on the Hubble 
Space Telescope. Finally, we summarise our findings in §5. 

We provide on-line tables including cooling and photo- 
heating efficiencies per ion, ionisation fractions in equilib- 
rium and non-equilibrium cases, and the compiled atomic 



2 METHOD AND TESTS 

Our method explicitly follows the ionisation states of all 
11 elements that contribute significantly to the cooling ef- 
ficiencies at temperatures T ^ 10 4 K (i.e. H, He, C, N, 
O, Ne, Si, Mg, S, Ca, & Fe; W09) as well as the electron 
density. These are all the species tabulated by W09 and 
used in the Overwhelmin gly Large Simulations (OWLS) 
project ([Schave et al.ll2Q10h. Thev account f or nearly all the 
radiative cooling (W09; Ber tone et al.ll2013h . Our method is 
meant to be complimentary to W09, so that we can turn off 
following the non-equilibrium ionisation states of any metal 
species and instead use the corresponding equilibrium ta- 
bles. For example, iron involves following 27 ionic species, 
which may exacerbate the memory requirements when run- 
ning a large simulation. The code can be used to follow only 
the metal-line coolants that dominate at a given temper- 
ature (C, O, Ne, & Fe), or it can be used to follow the 
ionisation states of species frequently detected in QAL ob- 
servations (e.g. C, O, Si, & Mg). Species including S and 
Ca, which have many ions but are not very important for 
cooling and rarely seen in QAL spectra may be less useful 
to follow, so they can be turned off. We also build in options 
to turn off Auger ionisation and charge transfer reactions, 
because these calculations can significantly increase compu- 
tation time in exchange for mostly minor changes in the 
ionisation balance. 

Having computed the ionisation states for each species, 
we can use ion-by-ion lookup tables as functions of temper- 
ature and redshift for the cooling and photo-heating effi- 
ciencies. Because we explicitly track ionisation states, the 
lookup tables are much more compact than those of W09, 
who had to track net cooling efficiencies as a function of den- 
sity, temperature, redshift, and helium fraction. Finally, we 
consider criteria to cycle and sub-cycle the time integration 
of ionisation and cooling. 

Throughout, we will assume the gas to be optically thin 
to both the EGB and the cooling radiation. Assuming the as 
to be optically thin is generally a good approximation f or gas 
with densities nu < 10 -3 cm -3 ([Rahmati et al.l 120120 , but 
we note that for higher densities self-shielding may reduce 
the effects of photo-ionisation. We will ignore sources of ion- 
ising radiation other than the EGB, including self-radiation 
from shocks and local sources, which is a conservative as- 
sumption when investigating the potential effects of photo- 
ionisation. We will also assume that the electron and ion 
temperatures are equal. Although the electron temperature 
may temporarily lag that of the ions in shock-heated gas, 
supernova remnants suggest that plasm a waves quickly equi- 
librate the different temperatures (e.g. iBertone et aU F2008. 
and references therein). 

We first discuss the method and tests of following the 
ionisation of the gas, then the cooling and photo-heating 
method and tests, and finally how we cycle and subcycle 
during timesteps. 
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2.1 Ionisation 

We consider the processes of radiative and di-electric recom- 
bination, collisional ionisation, photo- ionisation by an EGB, 
Auger ionisation, and charge transfer in our calculations of 
the ionisation balance. For an ionisation state i of element 
x, the time dependent evolution (omitting Auger ionisation 
and charge transfer for simplicity) of the number density, n, 
of each ion species, Xi, is given by 



dn Xi 



dt + 

- n Xi ((a Xi + p Xi )n e + T Xi 



+ n Xz _ 1 (/3 Xi _ 1 n e + 



(1) 



where n e is the free electron density (cm -3 ), a Xi is the total 
recombination rate coefficient (radiative plus dielectronic, 
cm 3 s -1 ), f3 Xi is the collisional ionisation rate coefficient 
(cm 3 s -1 ), and T Xi is the photo- ionisation rate (s -1 ) for 
the given ionisation state, where 



f 



4ttJ u 
hv 



. (y)dv, 



(2) 



where v is frequency, z/o,^ is the ionisation frequency, J v is 
the EGB radiation field (erg s _1 cm -2 Hz -1 sr -1 ), a Xi (y) is 
photo- ionisation cross-section, and h is the Planck constant. 

We use atomic data simila r to that used by CL OUDY 
ver. 10.0(0 (last described in iFerland et all 1 19981 ) where 
reasonably possible. For radiative and diel e ctroni c re- 
combination coefficients, we use the iBadnelll {2006) and 
iBadnell et al.l (|2003l ) fits, respectively, if available for 
our s pecies. Otherwise, we use the IVerner &; Ferlandl 
(1996) fits for radiative recombination if available, and 



finally recombination da ta comp ile d by D. Vernej 2 ] from 
lAldrovandi fe Peauignotl (Il973h: IShull fc van Steenberg 



(1982): lArnaud & Rothenflug (1985). We assume case A 
recombin ation. Fits to col lisional ionisation rate coefficients 
are from IVoronovl (|l99/t ). We tabulate these coefficients 
between log[T (K)]= 2.0-9.0 at 0.04 dex intervals resulting 
in 176 temperature bins, which is the same spacing as used 
by W09. 

We perform the integral in Equation [2] u sing the same 
(TxAv) as used in CLOU DY ver. 10.00 (|Verner fe Yakovlevl 



Il995l : IVerner et aj]ll996h for the appropr iate EGB. In this 



work, we use the lHaardt fe Madaul (|200lL hereafter HM01) 
EGB as our fiducial background since this results in good 
agreeme nt with the observed column den sity distribution 
of H I ([Dave et al.l l2010l; lAltav et al.l l201lh and metal- line 
ratios (jAguirre et al.l l2QQ8h . This EGB includes contribu- 
tions from quasars and galaxies and is tabulated at 50 red- 
shifts between z — — 9 sampling at constant log(l+z) 
intervals, resulting in a lookup table of r Xi s as a function 
of redshift corresponding to each ion containing electrons. 
We also tabulate photo-ion isat io n cross-sections for the 
iFaucher-Giguere et alJ (2009) and lHaardt fe M adau (2012) 
and make these available on our website. 

Auger ionisation involves the ejection of multiple elec- 
trons following the ionisation of an inner shell electron. The 
Auger reaction is stimulated by an energetic photon, 7, such 



that when element A, initially ionised to charge +n, inter- 
acts with 7, m electrons are released. The rate equation is 
then A +n + 7 A +n+m + m e~, where m > 1. We use 
the electron vacancy di stribution probabilities calculated by 
iKaastra fe Mewd (|l993h . and re-tabulated for CLOUDY as 
probabilities of ionising multiple numbers of electrons, all 
the way up to 10 electrons for the case of Fe I to Fe XI. Auger 
ionisation is significant for photo-ionised oxygen species, be- 
cause 2 electron losses are common for inner shell ionisations 
of O I through O V. We include the option to turn Auger ion- 
isation off, because it slows the calculation, or to lower m 
to values as low as two, if one does not want to follow the 
m > 2 reactions that are mostly limited to heavier species 
like Fe. 

Charge transfer reactions encompass two-body reac- 
tions where neutral H or He exchange an electron with a 
heavier ion species leading to recombination of the heavier 
ion. For hydrogen reacting with element A ionised to charge 
+n, the reaction is H° + A +n -+ H + + A +n ~\ The most 
significant reaction is the recombination of O II to O I owing 
to H 1, which can enhance O 1 levels. The reverse reactions 
of singly ionised H and He with a heavier species ionises the 
heavier species such that for hydrogen, H + + A +n_1 — ► H° 
+ A +n . CLOUDY uses r ates from various references (e.g. 
iKingdon fe Ferlandl 1 1996h to calculate the four varieties of 
charge transfer reactions: recombination and ionisation with 
H and He. We have outputted and tabulated these rates for 
our 176 temperatures from CLOUDY for all relevant metal 
ions with H and He, as well as charge transfer between H 
and He. 

We use the Sundials CV0D^| solver to integrate the ion- 
isation states over a timestep dt to solve for each xi. We use 
the backward difference formula (BDF) method and New- 
ton iteration in CV0DE, given the set of rate equations. For 
an atomic species with atomic number Z, there are Z + 1 
versions of Equation [1] to solve for the Z + 1 ions that exist; 
hence the number of equations equals the number of un- 
knowns. In practice we follow all Z + 1 ions, even though 
an independent constraint is available, namely that the sum 
of all ion fractions equals one, which allows one to follow 
only Z states. We calculate ^2 f x . for each element, where 

i 

f Xi s are ion fractions, as an independent check to see how 
well the BDF method is tracking the states, which depends 
on the accuracy that one chooses for the calculations using 
the relative and absolute tolerance parameters described in 
the CV0DE manual. Briefly, for the calculations in this pa- 
per we set absolute tolerance vector to 10 ~ 5 x f x . for metal 
species and 10 ~~ 7 x f Xi for hydrogen and helium. The rel- 
ative tolerances are 10 ~ 5 and are unlikely to play a role, 
since the absolute tolerance vectors are equal or tighter, plus 
they act as relative tolerances since they are multiplied by 
the ionisation fraction. However, for the version integrated 
into Gadget-3, we will choose looser tolerances according 
to the computational resources. If ^2 f Xi strays more than 

i 

1% from unity, we renormalize the ionisation fractions. How- 
ever, such deviations occur only in extreme situations and 
are especially rare given our low tolerances. Deviations from 
unity indicate that the timestep given to the CV0DE solver 



1 http://www.nublado.org/ 

2 http:/ /www. pa. uky.edu/~verner/rec. html 



https: / / computation, llnl.gov / case / sundials / main, html 
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is too large, which mostly affects the species with the most 
ions (e.g. Fe). 

We combine all atomic species into one CVODE calcula- 
tion, which includes Z + 1 ions for each atomic species, as 
well as n e , which we also follow in the solver. Including all 
eleven atomic species, we follow 134 rate equations. We in- 
clude the option to turn off the rate equation for n e if Auger 
ionisation and charge transfer are also turned off, because 
n e depends on as many as 133 ions and ruins the symmetry 
of the Jacobian matrix (which Auger ionisation and charge 
transfer also break), causing the calculation to scale non- 
linear ly with the number of rate equations. As most of the 
IGM is highly ionised and n e changes very slowly outside of 
the Hi and Hell reionisation epochs, we build in the option 
to update n e only after each CVODE iteration. For the cal- 
culations in this paper we follow everything self-consistently. 
Similarly, we leave Auger ionisation and charge transfer on 
in our calculations, but we generally recommend turning 
these off in a hydro simulation owing to longer computation 
times from non- symmetrical entries in the Jacobian. 

2.1.1 Ionisation test cases 

We demonstrate the capabilities of the ionisation solver by 
considering the evolution of initially neutral, solar abun- 
dance gas at a large range of densities and temperatures cor- 
responding to IGM/CGM regions. Figure [T] shows the ion- 
isation state of oxygen in density-temperature phase space 
at five distinct times (t = 0.1 Myr, 1 Myr, 10 Myr, 100 
Myr, & 1 Gyr) after the neutral initial conditions, and also 
in equilibrium, which is nearly equivalent to the case af- 
ter 10 Gyr0 In this example, we assume the z = 1 HM01 
EGB and neglect redshift evolution. We apply no cooling 
or photo-heating in this example, allowing us to focus on 
the ionisation behaviour by leaving the density and tem- 
perature constant. At t = (not shown), the entire phase 
space would be dark blue, corresponding to O I, while af- 
ter 1 Myr, the highest densities and temperatures (regions 
that would typically be defined as the intracluster medium 
(ICM), n H > 10" 3 cm" 3 , T > 10 7 K) have achieved com- 
plete ionisation to O ix, corresponding to dark red. 

The evolution of the oxygen ionisation state over the 
next Gyr demonstrates how significant ionisation timescales 
lead to non-equilibrium ionisation behaviour. The ionisation 
timescale for a given ion i is 

The collisional ionisation timescale, Tcoii,^ = VCA^e), 
depends inversely on density; hence the highest densities 
achieve equilibrium in less than a Myr. This includes gas 
collisionally ionised to O IX, which depends on eight ionisa- 
tion timescales, r ion q 1 through r ion Q vill' an °^ wn i cn are 
<C 1 Myr at the highest densities shown in Figure U For each 
1 dex increase in t, the density achieving collisional ionisa- 
tion equilibrium (CIE), defined as where collisional ionisa- 
tion exactly balances recombination, decreases by 1 dex. In 

4 For movies showing time evolution of the ionisation 
states, and also for other metal species, please visit 
|http://noneq.strw. le idenuniv.nl 



CIE, the ionisation state depends only on T, which is why 
there is no density dependence at the highest densities where 
photo- ionisation by the uniform HM01 EGB is unimportant 
(f3 Xi n e > T Xi ). In CIE the O VI fraction peaks at T « 10 5 5 
K, corresponding to the horizontal yellow band indicating 
this state on the right-hand side of each panel in Figure [T] 

The photo- ionisation timescales, Tphot,^ = l/T Xi , are 
independent of density for a uniform EGB. At low densities 
where Tphot,^ <C Tcoii,^ , the ionisation state depends on the 
series of r p hotS, which increase with ionisation state. After 
100 Myr the dominant ionisation state at low density (tih < 
10 -6 cm -3 ) is Ovi and after 1 Gyr it is Oviii. The lowest 
densities, corresponding to the centres of voids at z = 1 
(nu < 10~ 7 cm -3 ), take nearly 10 Gyr (assuming the z = 1 
HM01 field) to achieve their equilibrium ionisation states, 
which are dominated by O ix for all temperatures. 

The dashed vertical white line, indicates the density for 
which the dynamical timescale, 

Tdyn = 1/y/Gj, (4) 

equals t, where G is the gravitational constant and p is 
the gas mass density. Td yn is long compared to t to the 
left of this line. The line is beyond the right boundary of 
the ^10 Myr panels, meaning that the gas exhibiting non- 
equilibrium ionisation states will, at these times, not be dis- 
rupted by dynamical processes. An astrophysically relevant 
case may be the injection of enriched, low-ionisation wind 
material into the IGM, where an "ionisation la g" leaves the 
gas un der- ionised relative to equilibrium values. ICen &; Fangl 
(2006) found enriched gas was under-ionised in their cosmo- 
logical simulation that explicitly followed the oxygen states 
O V — O ix, which was reflected in fewer predicted O VII and 
Ovill absorption line detections in their mock AGN X-ray 
spectra probing the IGM. 

We provide a consistency check for the predicted ioni- 
sation fractions of our code by comparing to CLOUDY ver. 
10.00 calculations in Figure [2] for gas in equilibrium, which 
we obtain for our code by integrating for at least 100 x any 
recombination or ionisation timescale. 

The top panel shows the equilibrium ionisation fractions 
for Hi, Civ, Ovi, and Nevin as a function of temperature 
for a density where collisional ionisation dominates, nu = 
0.1 cm -3 . The dotted lines are from CLOUDY-generated 
ionisation lookup tables at z = 1 assuming equilibrium. Our 
code produces the solid lines that reproduce the CLOUDY 
CIE behaviour of the four ions, which depend on the balance 
of recombination and collisional ionisation. 

The bottom two panels of Figure [2] demonstrate the im- 
pact of photo-ionisation as a function of density for T = 10 4 
K (middle) and 10 5 ' 5 K (bottom). Photo-ionisation domi- 
nates in the T = 10 4 K case and the two codes agree well. For 
Hi, nm oc nu because recombination (oc n e ) balances photo- 
ionisation (independent of n e ) for this two- state species in 
equilibrium. In contrast, the lithium-like metal ions exhibit 
a peak ionisation fraction corresponding to a preferred den- 
sity. As the atomic number increases, the energy required 
for ionisation increases, and the number of photons capable 
of ionising decreases for realistic spectra. Hence, the density 
at which the equilibrium ion fraction peaks declines with 
increasing atomic weight. 

The T = 10 5 ' 5 K case shows that our code repro- 
duces the transition from photo-ionisation to collisional 
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Figure 1. The evolution of non-equilibrium oxygen ionisation fractions starting from a neutral state at t = for a range of densities 
and temperatures using the z = 1 HM01 radiation field. The colour indicates the average ionisation fraction, from O I (dark blue) to O ix 
(dark red). The first five panels correspond to the five distinct times indicated and the bottom right panel corresponds to the equilibrium 
case, which is reached by 10 Gyr for all densities. The vertical, dashed white line indicates where the dynamical time equals t, which only 
corresponds to densities in the displayed range at 100 Myr and 1 Gyr. This implies the area to the left of this line (the entire phase spaces 
at ^ 10 Myr) would not be dynamically disrupted. Collisional ionisation timescales are shorter at high densities, achieving equilibrium 
values earlier. Photo-ionisation timescales are independent of density and equilibrium has not yet been achieved by 1 Gyr at the lowest 
densities where photo-ionisation dominates. For movies of the evolution, please visit http://noneq.strw.leidenuniv.nl 



ionisation at higher densities for all species. We choose 
this temperature specifically because O VI achieves its CIE 
peak at this temperature, /hi asymptotes to 10 -6 in 
the collisional ionisation regime for this temperature. For 
further comparison tests between the equilibrium abun- 
dances predicted by our code and CLOUDY, we forward 
the reader to the Appendix and additionally our website- 
|http: / / noneq.strw.leidenuniv.nl , We also demonstrate the 
effects of Auger ionisation and charge transfer in the Ap- 
pendix. The examples in this Section have both of these 
processes turned on. 

While these comparisons demonstrate that our reaction 
network contains all the important processes and that our 
atomic data is close to that used by CLOUDY, they do 
not test the accuracy of time integration. We will test the 
time-dependent non-equilibrium behaviour in $3] when we 
compare the results of our method to GS07. 

2.2 Radiative cooling and photo-heating 

The net volumetric cooling rate (A ne t) is the difference be- 
tween cooling and photo-heating rates. W09 tabulated net 



cooling tables under the assumption of ionisation equilib- 
rium as functions of temperature, density, redshift, and he- 
lium fraction for 9 metals and primordial species separately. 
Here we make use of ion-by-ion cooling and photo-heating 
tables that are more compact and flexible, which can be 
combined into net cooling functions for both equilibrium 
and non-equilibrium calculations. Our ion-by-ion net cool- 
ing method uses the ion fraction calculated using the CVODE 
solver. We calculate net cooling as the difference between 
cooling rates (A^.n x .n e ) and photo- heating rates (e Xi n Xi ), 
where the resulting net cooling function for an ion x% is 

A ne t,x, (T, z, n x% , rte) = A.' x . (T)n x% n e - e Xi (z)n x% (5) 

The units of A' x . , ion-by-ion cooling efficiencies, are erg cm 3 
s _1 where we divide the volumetric cooling rate for an ion 
by n Xi n e to make A' independent of density (at least in the 
regime where electron-ion collisions dominate the cooling 
rate). Photo- heating is calculated in an analogous manner 
as photo-ionisation. The photo-heating efficiency for species 
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(6) 



and hence has unit ergs s _1 . We calculate e Xi for our 133 
Xi species in units of erg s _1 for every ion given the HM01 
EGB at the 50 redshifts between z = and 9, tabulating it 
in a table analogous to that of the photo-ionisation cross- 
sections. 

We have used our own method to calculate CLOUDY 
cooling efficiencies i on- by-ion (A^,). which was developed 
independently from I Gnat &; Ferlandl (|2012l ) , who also tab- 
ulate ion-by-ion cooling functions using the same version 
(ver. 10.00) of CLOUDY. We run CLOUDY to calculate in- 
dividual ion cooling efficiencies at 176 temperature values 
between log[T (K)]= 2.0 and 9.0. The end result is a set 
of cooling tables that include the effects of radiative and 
dielectronic recombination, collisional ionisation and excita- 
tion, and Bremsstrahlung. We have compa r ed ou r ion cool- 
ing functions to those of iGnat &; Ferlandl (|2012h and find 
that our efficiencies are indistinguishable from theirs 

In detail, we tabulate ion-by-ion cooling functions by 
running CLOUDY with the metal ion fraction we are calcu- 
lating set to unity. We assume solar abundances as was done 
in W09. We force hydrogen and helium to be fully ionised, 
resulting in n e = 1.2riH from these species alone. We then 
run CLOUDY again with the metal being calculated turned 
off, and difference the two sets of cooling efficiencies to ob- 
tain A' x . , which we tabulate for every metal ion in units of 
erg cm 3 s _1 , because we divide out n Xi n e from the volumet- 
ric cooling rates to maintain density invariance. In practice 
the difference in n e with the metal ion turned on and off 
is very small, since metal species do not contribute many 
electrons compared with H and He unless the metallicity is 
much higher than solar. The calculation for primordial ions 
is similar, except that the electron density will not default to 
n e = 1.2nu in the cases of Hi, Hei, and Hell, because a sig- 
nificant fraction of the electrons can become bound. There- 
fore, we force n e in our CLOUDY scripts to add in other- 
wise bound electrons for these primordial species, resulting 
in n e = 1.2 cm -3 . Our method differs from lGnat &; Ferlandl 
(|2012h . who essentially make a single ion plasma to calcu- 
late cooling efficiencies. A strategy similar to ours was used 
to create the W09 tables except that W09 did not force ion 
states since its cooling was calculated assuming ionisation 
equilibrium. 

The total summed net radiative cooling rate per unit 
volume of all species is 



A n 



+ Ago 



(7) 



Figure 2. Comparisons of equilibrium ionisation fractions for H I, 
Civ, Ovi, and Nevm predicted by our code and by CLOUDY 
using the z = 1 HM01 radiation field. The top panel shows ion 
fractions as a function of temperature at the constant density 
riH = 0.1 cm -3 , for which collisional ionisation dominates in most 
species. The bottom two panels show ion fractions as a function of 
density at T = 10 4 and 10 5-5 K, respectively. Photo-ionisation 
dominates the middle panel while the bottom panel tests the tran- 
sition between the photo- and collisionally ionised regimes. The 
agreement with CLOUDY is excellent, although there is devia- 
tion at some very low metal ion fractions where CLOUDY does 
not accurately track these very low abundances. 



where inverse Compton cooling off the cosmic microwave 
background contribute Acomp = 5.64 x 10 -36 erg K _1 
-l (r _ Tcmb.oKI + z) V , using Tcmb,o = 2.728 K 
(|Osterbrock fe Ferlandl 1 20061 ). 

The temperature evolution is calculated using 



dt ~ (3/2 + s)n/c ji dt ' 



(8) 



5 Further comparison figures are available on our website. 
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where n, the total number density of the gas, is p/(iimn) 
and /x, the mean particle mass, is defined as 



(9) 



nrriH 



where k is the Boltzmann constant, and s is for isochoric 
(i.e. constant density) cooling and 1 for isobaric (i.e. con- 
stant pressure) cooling. The isobaric case cools more slowly 
owing to performing "PdV" work, and we refer the reader 
to §2.2 and §3 of GS07 for a more complete explanation 
of isochoric and isobaric cooling. We mainly focus on iso- 
choric examples, because this is what is applicable to inte- 
grating cooling during a hydrodynamic timestep in a simula- 
tion. Lastly, the second term in Equation [8] accounts for the 
change in internal energy owing to the change in the number 
of particles sharing the internal energy in equipartition. 
We define the cooling timescale as 



''"cool — 



(3/ 2 + s)n/cT 

Anet 



(10) 



and provide this timescale with data included on our web- 
site. 

We will continue to update our ion-by-ion cooling data 
at our website http://noneq.strw.leidenuniv.nl We note the 
recent improveme nts in atomic data used to calculate ra- 
diative cooling by Lvkins et al.l (|2013h . Our website will be 
updated to provide new cooling and ionisation coefficients 
from the most recently released official version of CLOUDY. 



2.2.1 Cooling test cases 

The top panel of Figure [3] shows the cooling functions (A^.) 
for individual oxygen ion s and compares them to those 
of iGnat fe Ferlandi (|2012h . We check that the two meth- 
ods agree after scaling to CLOUDY default solar abun- 
dances The total volumetric cooling rate due to element 
x is ^2 A^.TixUefxi- Our CIE oxygen ionisation fractions 

i 

(solid lines) are compared to the CLOUDY fractions (dotted 
lines) in the second panel, demonstrating excellent agree- 
ment. The third panel shows the summed oxygen CIE cool- 
ing curve (thick black line) multiplied by n e /nu to plot 
cooling efficiency, combined with the contribution of each 
ion species (coloured lines). We compare this to the anal- 
ogous CLOUDY calculation (grey solid line), which agrees 
almost perfectly above T = 10 4 K and is virtually indis- 
tinguishable from our summed cooling curve. The cooling 
rate drops precipitously below T ~ 10 4 K, because in CIE 
n e declines steeply as hydrogen recombines. However, the 
CLOUDY calculation diverges and has much higher cooling 
owing to collisions other than electron-ion collisions, in this 
case H 1-0 1 collisions. Our ion-by-ion cooling method fails in 
this regime for CIE, but we can still safely apply our method 
to photo-ionised cases below T — 10 4 K where electron-ion 
collisions dominate cooling. 



These abundances, also used in W09, are nn e /nn = 0.1, 



10" 



n N /n H 



nc/nn 

%e/wH = 10~ 4 - 00 , n Mg /nu = 
ns/nu = 10 -4 - 74 , nca/nu = 10" 



10 



-4.07 



no/nu = 10 
: 10- 4 - 46 , nsi/nn = 10 
5 64 , and n Fe /n H = 10 
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Figure 3. The top panel shows the cooling efficiencies of in- 
dividual oxygen ions at solar abundance calculated by us using 
CLOUDY from T = 1 2 - 10 9 K (solid lines ) compared to ef- 
ficiencies calculated by IGnat fc Ferlandi (120121 ) using CLOUDY 
(dotted lines). The second panel compares the collisional ionisa- 
tion equilibrium (CIE) oxygen fractions we calculate using our 
reaction network (solid lines) to those predicted by CLOUDY 
(dotted lines). The integrated oxygen CIE cooling efficiency, plot- 
ted as A^ .Tie/riH (solid black line) in the third panel, is the sum 
of ion fractions multiplied by the ion cooling efficiencies (coloured 
lines), and agrees well with the CLOUDY case (grey line), except 
below T = 10 4 K where our method loses accuracy due to the 
presence of cooling channels other than electron-ion collisions. 
The photo-ionisation equilibrium (PIE) oxygen ionisation frac- 
tions and summed cooling efficiency appear in the bottom two 
panels, assuming nn = 10 -4 cm -3 and the HM01 EGB at z = 1. 
Cooling, photo-heating (e Xi /nn), and net cooling (A ne t/^H 2 ) 
indicated in the bottom panel by dashed, dotted, and solid lines 
respectively for our code (black) and CLOUDY (grey). The equi- 
librium temperature of the solar abundance gas is indicated by 
the vertical dashed line. 
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In the presence of a photo-ionisation background three 
things happen. First, ionisation removes electrons meaning 
higher ions exist at lower temperatures. This alone reduces 
cooling efficiencies in most cases, since higher ions usually 
have lower A' x . at a given temperature, because fewer ra- 
diative transitions exist; however there are exceptions. The 
fourth panel in Figure [3] shows the photo- ionization equilib- 
rium (PIE) ion fractions of oxygen for our code (solid lines) 
and the CLOUDY calculation (dotted lines). Here we use a 
density n H = 10~ 4 cm" 3 and the HM01 EGB at z = 1. For 
states T > 10 5 K the agreement is excellent while there are 
some small, but noticeable deviations at lower temperatures 
for O in, O IV, and O V. The PIE cooling efficiencies in the 
bottom panel show very good agreement. Our investigation 
in AppendixlAl shows that low ions often agree less well with 
CLOUDY calculations than high ions. 

Second, n e does not decline at T < 10 4 K, because 
photo-ionisation keeps most hydrogen ionised. Hence, for 
T < 10 4 K the summed cooling of oxygen in the bottom 
panel (dashed thick black line) is actually much higher for 
PIE than for CIE owing to the much greater electron density. 
We provide tables down to T = 10 2 K, because electron-ion 
collisions dominate the cooling in this regime given a realis- 
tic background. 

Third, the total volumetric photo-heating of oxygen, 
^2^xin x f Xi , represented by the dotted thick black line in 

i 

the bottom panel, adds energy to the gas. We do not show 
the contributions of individual species to the photo-heating 
as we do for the cooling. We plot net cooling, which is de- 
fined here as cooling minus heating, as the solid thick black 
line. For T < 10 4 K even the net cooling due to oxygen is 
greater in PIE than in CIE! The point where cooling equals 
heating is 10 2-50 K for oxygen. However, the thermal equilib- 
rium point for this entire solar abundance gas is at 10 4 ' 23 K, 
mainly owing to hydrogen, and is indicated by the dashed 
vertical line in the bottom two panels. We compare to the 
CLOUDY values for cooling, heating, and net cooling with 
thick grey lines, finding very good agreement. 



2.3 Cycling and sub-cycling 

Given a timestep, which we will call the cycle timestep, dt, 
which may be the Courant timestep in a hydrodynamic code, 
we determine if the cooling rate needs to be updated on a 
shorter, sub-cycle timestep by determining if Su/u > £ CO oi 
where u is the internal energy density, Su = A ne t x dt is 
the change in this energy, and £ CO oi is the fractional cooling 
tolerance. We then update the ionisation state explicitly on 
either the cycle timestep if Su/u ^ £ C ooi, or on the sub-cycle 
timestep defined by u/A ne t x £ C ooi- For the calculations pre- 
sented here, we use £ CO oi = 0.01, but we find that running 
with 0.05 gives acceptable accuracy for much less computa- 
tion time. 

If we need to sub-cycle, we iterate cooling and ionisa- 
tion, where the latter is integrated using CVODE as dis- 
cussed in fc|2.1l explicitly until either i) the cycle timestep 
is reached, or ii) the equilibrium temperature is achieved 
by noting that the sign of A ne t reverses. If a sign reversal 
occurs, we backtrack to the equilibrium temperature where 
Anet = 0, since further subcycling is unnecessary and compu- 
tationally expensive. This is mainly a cosmetic choice as the 



internal energy will wiggle around the equilibrium tempera- 
ture with a tolerance according to £ CO oi between cycles if the 
particle remains at the equilibrium temperature. The ionisa- 
tion state will also wiggle, because the ionisation timescales 
are typically very short; hence our extra timestep attempts 
to prevent this unsightly wiggling. 

Finally, it should be noted that the ionisation sub-cycle 
timestep is determined by cooling. The reverse, i.e. the cool- 
ing timestep being sub-cycled according to changes in the 
ionisation state, is not possible in our explicit method de- 
scribed thus far. In most cases applicable to diffuse gas, the 
cooling efficiency is not dramatically altered by a change in 
the ionisation state, but there are some notable exceptions. 
The most obvious example is reionisation, where the ionisa- 
tion timescale is short and the cooling characteristics change 
dramatically from neutral to ionised. Another case is that 
of a nearby AGN suddenly turning on or off, changing the 
ionisation state on a timescale much shorter than the cool- 
ing time. We will explore the latter case in Oppenheimer & 
Schaye (in preparation), where we will also discuss how we 
address the sub-cycling. For now, we state that our current 
method is meant to follow cooling in the presence of a slowly 
changing EGB. 



3 COLLISIONALLY IONISED 

NON-EQUILIBRIUM METAL-ENRICHED 
GASES 

Before considering non-equilibrium ionisation in the pres- 
ence of ionising radiation in we will first consider only col- 
lisional processes and compare to previous results from GS07 
(who used rates from CLOUDY ver. 06.02) and W09 (who 
used CLOUDY ver. 07.02) . Figure H shows O VI (left panel) 
and O in (right panel) ion fractions as a function of temper- 
ature for cooling gas initially at T = 10 7 K for the CIE and 
collisional non-equilibrium (CINe) cases at Z = 10" 3 ,0.1, 1, 
and 2 Z©. Our calculations (solid lines) are compared to 
GS07 (dotted lines), who show the same two ions in their 
Figure 3. We have scaled our metal abundances to the abun- 
dances in their Table 1 for direct comparison. These calcu- 
lations assume isochoric cooling, which is the default un- 
less noted otherwise. We start in ionisation equilibrium at 
T = 10 7 K, but the results are indistinguishable if the initial 
gas temperature is 5 x 10 6 K as in GS07, because cooling 
times above this temperature are long and non-equilibrium 
ionisation is unimportant. 

Figure 3] shows that the agreement with GS07 is good, 
and demonstrates the effects of non-equilibrium processes 
on ionisation fractions. Higher metallicities result in shorter 
cooling times and greater ionisation lags as recombination 
cannot catch up with cooling. Even for negligible metallic- 
ity, 10 -3 Z© (cyan lines), which is so low as to not impact 
cooling times relative to primordial, there is non-equilibrium 
behaviour for O in (right panel). Primordial species have 
non-equilibrium ionisation lags, and GS07 demonstrated 
that this behaviour is independent of metallicity below 
- 10~ 2 Z . The main culprit is the ionisation lag of re- 
combining hydrogen. Finally, it is worth noting that charge 
transfer is important for O ill. Without recombination of 
O in owing to charge exchange, O in would persist at tem- 
peratures much below 10 4 K. The slight mismatches between 
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Figure 4. A comparison of the predicted ion fractions of O VI (left panel) and O ill (right panel) as a function of temperature and 
computed assuming pure collisional ionisation (solid lines) with those of GS07 (dotted lines, see their Figure 3). The gas is assumed to 
cool isochorically, starting out in ionisation equilibrium at T = 10 7 K. The collisional ionisation equilibrium case (blue) shows excellent 
agreement between the two methods, as do the non-equilibrium cases (other colours corresponding to different metallicities). Even for 
low metallicity, Z = 10 -3 Z (cyan), non-equilibrium effects of primordial species cause O ill ionisation fractions to diverge from CIE. 
Metal abundances have been scaled to the values assumed by GS07. 



our method and GS07 are likely due to the older atomic data 
used in CLOUDY ver. 06.02 by GS07. 

Figure [5] shows the predicted isochoric cooling efficien- 
cies (solid lines), plotted as A ne t /(nune) for CIE (blue for 
10~ 3 Z©, magenta for Z ) and CINe (cyan for 10~ 3 Z , 
orange for Z©). We compare to W09 (dashed lines) for 
the CIE cases, and to GS07 (dotted lines) in all cases. We 
use W09 solar abundances here and throughout the rest of 
the paper and note that the GS07 abundances are differ- 
ent. The differences with respect to GS07 for the cooling 
efficiencies at T ~ 10 5 — 10 6 K is attributable to different 
neon abundances, which are nearly twice as high for GS07. 
The increased cooling at T > 10 6 K for GS07 is due to up- 
dated atomic data between CLOUDY ver. 06.02 and ver. 
10.00 (jGnat fc Ferlandl l2012h . These atomic updates were 
included in the CLOUDY ver. 07.02 used to generate the 
W09 tables, which show excellent agreement with our CIE 
cooling. 

In CIE, the 10 -3 Z curves have peaks below 10 5 
K corresponding to the collisional excitation of hydrogen 
and helium, and a drop off at higher temperatures until 
Bremsstrahlung dominates above 10 6 K. Solar abundances 
increase cooling efficiencies by more than an order of magni- 
tude between T - 2 x 10 5 and 2 x 10 6 K. Below T = 10 4 K, 
our method loses accuracy because cooling is no longer dom- 
inated by electron-ion collisions for this purely collisional 
example. The W09 Z cooling efficiency jumps up, because 



ion collisions with neutral hydrogen dominate the cooling 
here. However, the cooling efficiency is much lower below 
T — 10 4 K than the W09 magenta efficiency suggests, and 
the jump in the W09 rate is because we are dividing by n e 
(to compare with GS07), which is very small. Hence, the 
cooling rate A ne t is very small below 10 4 K both for W09 
and us. Nevertheless, care should be taken when our method 
is applied to situations where cooling is not dominated by 
electron-ion collisions. 

The CINe cases show the same characteristics as the 
corresponding GS07 efficiencies. In non-equilibrium the high 
ionisation fractions for primordial and metal species persist 
down to lower temperatures, where they are generally less 
efficient coolants. CINe cooling curves therefore are both 
smoother and lower than the corresponding CIE curves. This 
effect is also important for primordial species, which have 
very peaked CIE cooling efficiencies. 

Below T = 10 4 K, our method works for CINe cases 
as long as electron-ion interactions dominate the cooling, 
which is the case for solar metallicity. Vasiliev ( 2013) explore 
CINe cases including molecules, finding that for Z > 0.1 Z 
molecular cooling is negligible. Metals provide free electrons 
below T < 10 4 K, which are far more efficient coolants than 
molecules or ion-ion interactions. 

We compare the ionisation fractions for three collisional 
cases and solar metallicity in Figure [6] CIE Isochoric (solid 
lines), CINe Isochoric (dashed lines), and CINe Isobaric 
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Figure 5. Isochoric cooling efficiencies (plotted as A ne t / '(nnn e )) for collisional ionisation equilibrium (CIE) and collisional ionisation 
non-equilibrium (CINe) (solid lines) compared with W09 (dashed lines) and GS07 (dotted lines) for 10 -3 Zq and Zq. The agreement 
with W09 is excellent everywhere above 10 4 K. Comparison of our efficiencies to GS07 shows small deviations at higher temperatures 
partially owing to the different solar abundances used by GS07. We find the same non-equilibrium behaviour as GS07 of smoothed cooling 
curves as ion states persist down to lower temperatures where they generally are less efficient coolants. CINe curves begin in equilibrium 
at T = 10 7 K, but CIE and CINe curves are indistinguishable at T > 10 6 K. 



(dotted lines). For the isobaric case, we normalise the density 
to the isochoric cases at T = 10 6 K. A variety of observa- 
tionally obtainable metal ions are shown, starting with car- 
bon species on top, oxygen species in the second panel, and 
silicon species combined with the Li- like ions Nv, Nevin, 
and Mg x in the third panel. Our ion fractions compare well 
with GS07 (their Figure 2), exhibiting the behaviour ex- 
plained in their Figure 4. Briefly, He-like ions (C V, O vil) 
show the greatest ionisation lag, owing to having the longest 
recombination times. Isobaric cases do not have as great of 
an ionisation lag as isochoric cases, because of their longer 
cooling times at the same density due to extra "PdV" work. 
Our CIE ionisation fractions diverge from GS07 at T ^ 10 4 



K (cf. C II and Si n to their Figure 2) likely owing to ion- ion 
collisions, which are not included in our method, keeping the 
gas more ionised at T < 10 4 K. 

Finally, non-equilibrium effects are most important 
when cooling is most efficient, which is one reason we plot 
cooling times in the bottom panel. We use the specific case 
of 71h = 10 -4 cm -3 gas for comparison later with the photo- 
ionisation case in ^Hthat depends on the density in a non- 
trivial manner, but these cooling times scale inversely with 
7iH, which is why we plot t CO o1^h/(10 _4 cm -3 ). For CINe, 
cooling is less efficient above 1.3 x 10 4 but it is much more 
efficient at T < 10 4 K, because n e is much higher owing to 
H II persisting down to lower temperatures due to the recom- 
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Figure 6. Collisional-only ionisation fractions for the same cases explored in GS07: collisional ionisation equilibrium (CIE Isochoric, solid 
lines), collisional ionisation non-equilibrium at constant density (CINe Isochoric, dashed lines) and at constant pressure (CINe Isobaric, 
dotted lines). We show ionisation fractions for some of the most common species observed in diffuse gas in the top three panels, as well as 
cooling times in the bottom panel to indicate how long cooling gas lives at each state, assuming a density of nu = 10 -4 cm -3 and solar 
metallicity. The CINe Isobaric case matches the density of the CIE and CINe Isochoric cases at T = 10 6 K (i.e. nn = 10 -4 cm _3 (10 6 
K)/T). The CIE Isobaric case is not shown but would have the same ionisation fractions as the CIE Isochoric case while r coo i would be 
a factor 5/3 larger due to "PdV" work. 



bination lag, and additionally H n is a more efficient coolant 
than H I at these temperatures. The CINe Isobaric case has 
density nu — 10 -4 x (10 6 K/T) cm -3 , which means that 
for the same parcel of gas, the cooling will become faster at 
lower temperatures as its density rises. However, the cool- 
ing efficiency (erg cm 3 s _1 ) as plotted in Figure [5] is lower 
for the CINe Isobaric case compared to the CINe Isochoric 
below T — 10 4 K owing to quicker recombination to metal 
ions that are less efficient c oolants. This behaviour is also 
identified by Vasiliev (2013) (cf. their Figure 5 for Z©). 



4 PHOTO-IONISED NON-EQUILIBRIUM 
METAL-ENRICHED GASES 

We now explore non-equilibrium processes in photo-ionised 
gases, which is the main innovation in this paper. Our mo- 
tivation stems from the reality that metal-enriched diffuse 
gases cool in the presence of an ionising background, and 
that the cooling efficiencies, as well as the observational di- 
agnostics, are altered by the addition of photo-ionisation 
to the non-equilibrium processes described in the previous 
section. The density-independence of collisional processes is 
broken, because the photo-ionisation rate per unit volume 
due to a uniform field scales as n, while collisional processes 
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Figure 7. Ion fractions as a function of temperature for O VI (top left) and C IV (top right) for 6 cases at solar abundances: collisional 
ionisation equilibrium (CIE Isochoric, solid cyan), collisional non-equilibrium isochoric (CINe Isochoric, dashed cyan) and isobaric (CINe 
Isobaric, dotted cyan), photo-ionisation equilibrium (PIE Isochoric, solid red), photo-ionisation non-equilibrium isochoric (PINe Isochoric, 
dashed red) and isobaric (PINe Isobaric, dotted red). We also show the time evolution from T = 10 6 K at t = yr in the bottom four 
panels of temperature, electron density, O VI and C IV fraction. We assume equilibrium initial conditions of nn = 10 — 4 cm -3 and T = 10 6 
K using the z = 1 HM01 field in photo-ionised cases. All cases achieve T < 10 5 K temperatures in under 1 Gyr with cooling trajectories, 
electron densities, and ionisation fractions as a function of time differing significantly. 



scale as n 2 . Our fiducial case is a nu = 10 -4 cm -3 , solar 
abundance gas ionised by the HM01 EGB at z = 1. This 
corresponds to an overdensity, 5 ~ 60, and we will show 
that gas with this density cools in less than a Hubble time. 



4.1 Non-equilibrium behaviour 

We consider the three cases explored in the previous section 
(CIE Isochoric, CINe Isochoric, and CINe Isobaric) along 
with three new situations with the added z = 1 HM01 
EGB: PIE Isochoric, photo- ionised non-equilibrium (PINe) 
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Figure 8. Cooling efficiencies, plotted as A ne t/^H 2 , for the six cases displayed in Figure [Jlassuming solar abundances at nu = 10 -4 cm -3 
and the HM01 z = 1 EGB. The isobaric cases are set to nn = 10 -4 cm -3 at T = 10 6 K. W09 cooling curves are displayed as thin dashed 
black lines for CIE and PIE at z = 0.957, and agree very well with our equilibrium models. Photo-ionisation suppresses cooling by removing 
bound electrons, most importantly from Hi. Non-equilibrium effects at constant density suppress cooling owing to recombination lags 
where higher ions persist to lower temperatures where they are less efficient coolants. The isobaric cases follow the CINe Isochoric cases 
relatively closely, since the density dependence is divided out, although the PINe Isobaric case approaches an equilibrium temperature 
owing to photo-heating just below T = 10 4 K and just above tih = 10 — 2 cm -3 . Non-equilibrium effects on cooling are smaller in the 
presence of ionising radiation, and the effect of photo-ionisation is smaller for isobaric than for isochoric cooling. 



Isochoric and PINe Isobaric. These six cases are shown in 
Figure [3 The gas begins in all cases in ionisation equilibrium 
at T = 10 6 K and nu — 10 -4 cm -3 . The top two panels show 
O VI and C IV fractions as functions of temperature. Below 
this temperature, the cooling times become short enough to 
lead to non-equilibrium behaviour. 

Comparing PIE to CIE, for PIE the O VI and C IV ion 
fractions are much higher below 2 x 10 5 K and 6 x 10 4 K, 
respectively, while these fractions are lower around the CIE 
peak. All photo-ionised curves asymptote to the equilibrium 
temperature where A ne t = 0. Hence, if gas cools to the equi- 



librium point and stabilises (e.g. if Td yn ^> T CO oi), then the 
gas will maintain the ionisation fractions corresponding to 
the equilibrium temperature. 

Non-equilibrium effects in general allow ions to per- 
sist to lower temperatures due to significant recombination 
times. For the PINe cases this effect is reduced as photo- 
ionisation also enhances the degree of ionisation in the equi- 
librium case. The impact of the photo-ionising background 
is smallest for the isobaric case because the density increases 
when cooling at constant pressure, thus decreasing the ion- 
isation parameter. 
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The temperature evolution in the top of the lower panels 
of Figure [3 shows the unique thermal histories of gas in the 
six cases we explore. These behaviours are also reflected in 
the cooling curve plot in Figure [8] where we plot A ne t /nu 2 
to show net cooling efficiency for the six cases. We chose 
71h — 10 -4 cm -3 at z = 1 to show a case where photo- 
ionisation has a significant effect in gas that will cool in less 
than a Hubble time. 

Photo-ionisation suppresses cooling relative to the 
collisional ionisation case whether in equilibrium or 
non-equilibrium. Considering the isochoric cases, non- 
equilibrium effects retard cooling compared to the equilib- 
rium case resulting in the least efficient cooling curve of 
the six cases (Figure [8]) and the hottest gas after the peak 
of the cooling efficiency (T < 10 5 K) at late times (Fig- 
ure (T)) . The non-equilibrium suppression is greater for the 
collisional-only cases. Since the gas asymptotes to a highly 
ionised equilibrium state in the presence of ionising radia- 
tion, the non-equilibrium recombination lag affecting n e and 
riHi plays a smaller role and the cooling efficiencies of PIE 
and PINe Isochoric become nearly the same below T = 10 4 ' 5 
K (for riH = 10 -4 cm -3 at z = 1). Nonetheless, the reduced 
peak cooling rate of the photo-ionised cases delays the cool- 
ing of the gas by ~ 10 2 Myr in this example. The isobaric 
non-equilibrium curves also show differences, but the effect 
of photo-ionisation is smaller, because these cases achieve 
higher densities at temperatures below T = 10 6 K mak- 
ing photo-ionisation less important. The largest difference 
between the CINe and PINe isobaric cases is that the lat- 
ter asymptotes to an equilibrium point after just 800 Myr 
when Anet = 0, which is achieved at a lower temperature, 
T = 10 3 - 8 K, for the isobaric density, n H = 10" 1 ' 8 cm" 3 . 

Figure [8] also shows that we reproduce the W09 PIE 
cooling curve very well with our ion-by-ion cooling scheme 
(cf. black thin dash-dot line and red solid line). 

Figure [9] shows the dependence of the PINe Isochoric 
evolution on the density. We scaled the time axes by tih -1 , 
normalised to 10 -4 cm -3 , to distill the effects of photo- 
ionisation on the temperature and ionisation. The CINe Iso- 
choric case is shown in black for comparison. The effect of 
photo-ionisation is smallest at the highest densities. Lower- 
ing the density results in more delayed cooling even after re- 
moving the 71h -1 dependence of the cooling time, as well as 
higher equilibrium temperatures (i.e. A ne t = 0) because the 
photo-heating rate decreases less rapidly with the density 
than the cooling rate in this regime. The nu = 10 -5 cm -3 
case has the most extreme suppression of cooling shown, 
yielding a cooling time from 10 6 to 10 5 K of 19 Gyr, i.e. 
longer than a Hubble time. The EGB would therefore evolve, 
becoming weaker below z = 1, but we have turned off this 
evolution for simplicity. Finally, we note that Figure [9] also 
shows the effect of adjusting the strength of the HM01 field 
by a factor of 10 -4 cm _3 /riH at fixed density if we keep 
the time axis fixed in absolute terms. For example, the blue 
7i h = 10 -5 cm -3 curve is applicable to a nu = 10 -4 cm -3 
gas cooling in a 10 x stronger field and down to 10 5 K in 1.9 
Gyr. 

To summarise, photo-ionisation reduces cooling efficien- 
cies and increases the degree of ionisation for both equilib- 
rium and non-equilibrium cooling. Non-equilibrium effects 
on the evolution of the cooling rate and on the ionisation 
balance are smaller in the presence of ionising radiation. 



PINe Isochoric log[T] = 6, z=l, Z = Z Q 
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Figure 9. Evolution of non-equilibrium isochorically cooling solar 
abundance gas at various densities exposed to the z = 1 HM01 
radiation field (coloured lines) and assuming collisional ionisation 
only (solid black line). Ion fractions of Ovi and C IV are shown 
along with the net cooling efficiency and temperature. The time 
axis is scaled by wh and normalised to nu = 10 — 4 cm -3 to distill 
the effects of a uniform radiation field, which suppresses cooling 
more at lower densities. The PINe cases also illustrate the effect of 
turning up or down the HM01 field strength at tih = 10 -4 cm -3 , 
e.g. the blue curve represents a 10 x field increase and the red 
curve a 100 x decrease. 

The effect of photo-ionisation is smaller for isobaric than for 
isochoric cooling. 

4.2 Metal-line diagnostics 

QAL metal-line observations are used to diagnose both the 
physical state of the gas (e.g. its temperature, ionisation 
parameter, and chemical composition) and the dynamical 
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state of the gas (e.g. radiative cooling, conductive inter- 
face, turbulent mixing layer, shock front). Calculations usu- 
ally assume equilibrium ionisation or ignore the effect of 
photo-ionisation, which is why we emphasise the need to 
include both non-equilibrium and photo-ionisation in diag- 
nosing metal lines. 

The bottom two panels of Figure demonstrate the evo- 
lution of the O VI and C IV ion fractions for our six fiducial 
cases and the bottom panels of Figure [9] show the evolution 
for different densities assuming PINe Isochoric cooling. Note 
that for radiatively cooling gas, the probability of observing 
the ion ratios at the peak of the cooling curve is small com- 
pared to the probability of observing the gas at either higher 
temperatures where cooling is slower, or at lower tempera- 
tures when the gas approaches thermal equilibrium when 

Anet 0. 

In Figure fTQl we plot ion fractions for a variety of ob- 
served UV resonance transitions observed in QAL spectra 
for the PIE and PINe Isochoric cases for gas with nn — 
10 -4 cm -3 that is exposed to the z — 1 EGB, as well as the 
CINe Isochoric case. Allowing for the presence of ionising ra- 
diation opens up a new set of solutions with this figure show- 
ing only one density for a specific EGB. Entirely new fami- 
lies of solutions are possible for varying densities and back- 
grounds, which is why we offer the ionisation fraction tables 
for several densities for the HM01 background at several red- 
shifts on our website http: //noneq .strw. leidenuniv.nl , and 
inclu de the solutions also for the newer lHaardt &; Madaul 
(2012) background. Additionally, we provide cooling efficien- 
cies, Anet /mi 2 , an d the Tcooi as a function of temperature 
for various densities for the PIE Isochoric, PINe Isochoric, 
and PINe Isobaric cases. 

Figures El and fTQl show that assuming CINe rather 
than PINe can easily lead to order of magnitude errors in 
the predicted ion fractions. An example of an application for 
which the difference is much smaller is a QAL component 
absorber with A^NeVi n ^ N^y. This is very close to what 
has been observed bv I Tripp et al.l (| 20 111 ) who constrain the 
temperature of a radiatively cooling gas based primarily on 
these two species having similar columns (within 0.5 dex of 
each other) in various components of a single system. Since 
the solar number density abundances of nitrogen and neon 
are similar, A^Nevm ~ A^nv means that /Nevin ~ /nv- The 
rapid change in the N v and Ne viii fractions as a function of 
temperature visible in the third panel of Figure [TO] provides 
a powerful constraint for a single-phase radiatively cooling 
gas model of T ~ 10 5-5 K. By seeing where the respective 
green and red lines cross in this panel, we can see that photo- 
ionised models predict only slightly lower temperatures (T = 
10 5 ' 46 K for PINe Isochoric) than for the collisional-only case 
(T = 10 5-58 K). The ionisation correction also decreases for 
the PINe case as the ionisation fractions are 0.2 dex higher 
than for CINe meaning that taking photo-ionisation into 
account yields a lower gas mass. 

However, we provide a cautionary note about applying 
a single-phase radiatively cooling solution to this situation, 
because the probability of catching gas cooling at a tem- 
perature where iVNeViii ~ A^nv is low owing to r CO oi be- 
ing short (bottom panel of Figure 110 j) . Gas at these tem- 
peratures can exist, but it is a priori much more prob- 
able to observe cooling gas at either lower temperatures 
where r CO oi approaches infinity at the equilibrium tempera- 



ture, or at T > 10 6 K where r CO oi is longer. A multi-phase 
model comprising of radiatively cooling sub-clumps occupy- 
ing a temperature distribution weighted by t coo i (T) between 
T ~ 10 4 - > 10 6 K is more conceivable. However, this model 
would fail likely owing to the hotter sub-clumps having sig- 
nificant Mg x abso rption violat ing; the upper limit constraint 
observed bv lTripp et al.l ([201 ll ). A related possibility is that 
metal-enriched gas has a maximum of T~ 10 5 - 5 K accompa- 
nied by a much larger reservoir at lower temperatures, which 
would not violate the Mgx upper limit. However, we pro- 
vide an alternative model in Oppenheimer & Schaye (2013), 
which explains the observations by appealing to enhanced 
fossil photo-ionisation from a local AGN. 



4.3 Ionisation background variations 

Cooling efficiencies and QAL diagnostics both depend on 
the spectral shape and strength of the EGB. Thus far, we 
have used the HM01 EGB because it is wid ely used (e.g. 
ISchave et al.l l2010l : lOppenheimer et all 120121 ) and because 
it fit s observed metal ion ratios at z ~ 3 ([Aguirre et al.1 
2008) , as well as th e full Hi column density distribu- 
tion (jDave et alJl20ld : lAltav et al.ll201lh . However, the in- 
tensity and especially the shape of the EGB are rather 
uncertain as they depend on poorly constrained escape 
fractions, intrinsic spectral slopes, and extrapolations of 
luminosit y functions. Recently p ublished backgrounds in- 
clude the | Haardt fe Madaul (|20ll hereafter HM12) and the 
iFaucher-Giguere et al.1 (|2009l . hereafter FGllQ) model. We 
plot these backgrounds at the output nearest to z = 1 in the 
top panel of Figure [TT] to demonstrate the variations in the 
EGB, especially at the extreme UV (EUV) energies that cor- 
respond to the ionisation potentials of commonly observed 
Lithium-like ions indicated in the Figure. The assumed QSO 



EUV slope (>13.6 eV) ranges from zy -15 to zy -18 in the 
above listed models and can change the strength of the field 
by more than an order of magnitude at the ionisation energy 
of O vi (114 eV). Given that the observe d QSO EUV slope 
ranges between at least z/ -3 and z/°'° (jTelfer et al.l l2002i : 
IShull et aDl2012l ). we generate an EGB with a harder input 
QSO EUV slope, z/ -10 , using a post-2001, pre-2012 version 
of the CUBA package (F. Haardt, private communication). 
This background (H-1.0) is shown as the dotted blue line in 
Figure 1111 This EGB may be relevant for diffuse metals in 
regions occupied by harder than normal QSOs. 

Considering our fiducial case (tih = 10 ~ 4 cm -3 at z = 1 
and solar abundances), we show the PINe Isochoric cooling 
efficiencies for these four backgrounds in the middle panel 
of Figure fTTI The cooling time from 10 6 10 5 K, i.e. near 
the peak of the cooling curve, is 675 Myr for CIE and 773 
Myr in CINe. The addition of photo-ionisation from the 
three published EGBs increases this cooling time to between 
820 and 903 Myr. The more extreme H-1.0 EGB heavily 
suppresses metal coolants at these temperatures, increasing 
Tcooi to 1.47 Gyr. Hence, the inclusion of non-equilibrium ef- 
fects and photo-ionisation from a hard EGB can double the 
cooling timescale. Since the accretion of gas depends on its 



7 Using the Dec. 2011 update found 
https:/ /www. cfa.harvard.edu/~cgiguere/UVB. html 
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Figure 10. Ionisation fractions as a function of temperature for some of the most common species observed in isochoric, radiatively 
cooling diffuse gas along with cooling times (bottom panel) to indicate how long cooling gas lives at each state assuming wh = 10 -4 cm -3 . 
We assume an initial temperature of T = 10 6 5 K. The density is fixed at nu = 10 -4 cm -3 at all times. We assume solar abundances 
and, for the photo-ionised case (solid and dashed curves) the presence of the z = 1 HM01 EGB. In contrast with the analogous plot for 
collisional-only ionisation (Figure [6] and dotted curves in this Figure), a given ionisation species will often persist to lower temperatures 
due to photo-ionisation and stop at T = 10 41 K owing to photo-heating balancing cooling. For isobaric cooling (not shown here) the 
differences between PINe and CINe are smaller than for isochoric cooling. 



ability to cool, the process of galaxy formation can be sig- 
nificantly altered by photo-ionisation. This conclusion will 
only be strengthened if local sources of ionising radiation 
increase the intensity of the radiation field above that of the 
EGB. Such effects have been argued to be import ant (e.g. 
ISchavdl2006l : ICantalupollioTol : iRahmati et alll2013h . but we 
have conservatively ignored them here. 

While the most extreme EGB reduces cooling efficien- 
cies by a factor of two, the diagnostics of observed QAL 
can change much more as illustrated by the O VI fractions 
in the bottom panel of Figure 1111 Below the peak of the 



CIE cooling curve (i.e. at T < 10 5 3 K), Ovi fractions 
are much higher for PINe than for either CIE or CINe. 
The difference between the different PINe cases is smaller 
2x). We emphasise once again that, a priori, radia- 
tively cooling gas is most likely to be observed when r CO oi is 
longest, i.e. at the temperatures near the equilibrium point 
(T - 20, 000 - 30, 000 K). 
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Figure 11. Demonstration of how different ionising backgrounds 
(top panel) affect radiative cooling efficiencies (middle panel) and 
ion fractions (O VI shown in bottom panel). Different versions 
of the EGB at z = 0.957 show how changing the EUV slope 
can dramatically alter the ionising background at frequencies of 
commonly observed Lithium-like ion potentials. Non-equilibrium 
isochoric cooling curves for solar metallicity gas with density 
n H = 1CT 4 cm -3 beginning at T = 10 6 5 K at z = 1 (plotted 
as Anet/^H 2 ) show that harder radiation fields suppress cooling 
more by ionising metal coolants to a higher degree. For all back- 
grounds the cooling efficiency and O vi fraction are very different 
from the CIE and CINe cases (also shown in the bottom two 
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5 SUMMARY 

Radiative cooling is an essential ingredient for any model 
of the formation of galaxies and the evolution of the inter- 
galactic medium. Cooling efficiencies depend on the chem- 
ical composition and the ionisation balance. Knowledge of 
the ionisation balance is also critical for the interpretation 
of observations of quasar absorption lines. 

Although simulations nearly always assume ionisation 
equilibrium, this assumption will generally break down in 
gas that is cooling rapidly or that is exposed to a fluctuating 
radiation field. In particular, recombination lags may leave 
the gas over- ionised, thus reducing the cooling efficiency (e.g. 
lKafatod[l973l : IGnat Sternberdl2007h . 

Diffuse gas is exposed to the extra-galactic background 
radiation and possibly also to radiation from local sources. 
The presence of UV and X-ray radiation can strongly boost 
the ionisation of the gas and must therefore be taken 
into account when interpreting observations. The result- 
ing removal of bound electrons injects energy (i.e. photo- 
heating), but also suppresses the cooling efficie ncies of gas 
hot enough to cool through collisional excitation ( Efsta thioul 
Il992b IWiersma etHI l2009aL i.e. T > 10 4 K). Simulations 
have only recently begun to take th e se effects into acc ount 
for metals (e.g. (Schave et aD l2010l : iFord et al.1 l2012h , al- 
though still under the assumption of ionisation equilibrium. 

Clearly, a realistic treatment of ionisation and cooling 
can neither assume ionisation equilibrium nor ignore the ef- 
fects of the ionising background. In addition, the effects of 
non-equilibrium and photo-ionisation must be considered for 
both primordial and heavy elements. 

We have presented a method for calculating non- 
equilibrium ionisation and cooling of diffuse gas exposed to 
ionising radiation. We consider 133 ionisation states of the 
11 species that dominate the cooling efficiencies of diffuse 
gas (H, He, C, N, O, Ne, Mg, Si, S, Ca, & Fe). Cooling 
and photo-heating efficiencies are computed ion-by-ion us- 
ing tables generated with CLOUDY. After each time step, 
the ionisation balance is updated using a reaction network 
that includes radiative and di-electric recombination, col- 
lisional ionisation, photo-ionisation, Auger ionisation, and 
charge transfer. We make available the equilibrium and non- 
equilibrium calculations presented here in table format at 
http:/ / noneq.strw.leid enuniv.nl[ We also provide ionisation, 
recombination, cooling, and photo-heating coefficient input 
tables that can be integrated into hydrodynamic simula- 
tions, and explain our method to use these in such calcu- 
lations. 

We demonstrated that our method reproduces 
published collisional ionisation equilibrium and non- 
equilibrium ionisation f ractions and cooling efficiencies 
(|Gnat &; Sternberg! l2007h, and photo-io nised equilibrium 
ionisatio n fraction s (Ferla ncfet al.l [l998) and cooling effi- 
ciencies (Wiersma et a 1. 20093). Non-equilibrium ionisation 
and cooling in the presence of an extra-galactic background 
(EGB) is the main innovation presented here. 

Like IGnat &; Sternberg! (|2007r h we found that non- 
equilibrium cooling is confined to T < 10 6-7 K. At these 
temperatures, the recombination time becomes significant 
compared to the cooling time, resulting in a recombination 
lag and reduced cooling efficiencies. 

Photo-ionisation and non-equilibrium effects both tend 
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to boost the degree of ionisation and to reduce cooling effi- 
ciencies, sometimes by significant factors. The effect of the 
EGB is larger for lower densities (i.e. higher ionisation pa- 
rameters). Hence, photo- ionisation affects (equilibrium and 
non-equilibrium) cooling more under isochoric than under 
isobaric conditions. Non-equilibrium effects are smaller in 
the presence of an EGB and may therefore have been over- 
estimated by previous work. 

Observational diagnostics of diffuse, metal-enriched 
gases (e.g. metal absorption lines probed in quasar sight 
lines) are altered even more significantly than the cooling ef- 
ficiencies, both under equilibrium and non-equilibrium con- 
ditions. Hence, the inclusion of the EGB opens up an entire 
new range of ionisation solutions depending on gas density, 
metallicity, and the strength and shape of the ionising radi- 
ation field. 

We emphasised that not all ionisation solutions of a 
cooling gas should be treated with equal probability. For 
example, applying an ionisation solution to an observed ab- 
sorption system that corresponds to a temperature near the 
peak of the cooling curve (i.e. where cooling times are short- 
est) assumes all the gas traced by the absorber has the solu- 
tion that is a priori least probable. The addition of an EGB 
photo-heats the gas, usually leading to an equilibrium tem- 
perature where net cooling, radiative cooling minus photo- 
heating, is zero. Since the cooling time is effectively infinite 
at this point, this equilibrium solution with T ~ 10 4 K, has 
a high probability. 

Considering the dynamics of cooling is critical when 
applying a time-dependent solution to metal absorbers, 
which often require multi-phase models consisting of multi- 
ple parcels of gas at different temperatures. Hence, we pro- 
vide dynamical cooling times and efficiencies with all the 
ionisation fractions included on our website. Our fiducial 
case of solar abundance gas at nu = 10 ~ 4 cm -3 irradiated 
by a z = 1 EGB predicts much higher C IV and O VI fractions 
than pure collisional ionisation models. 

We include coefficient lookup tables extending from 
T = 10 2 — 10 9 K. Our method is applicable to cases when 
electron-ion collisions dominate the cooling, which includes 
nearly all intergalactic and circumgalactic gas. We demon- 
strate that our method works for metal-enriched gas at 
T < 10 4 K that is ionised by the EGB. However our method 
is not accurate for T ^ 10 4 K if H I is nearly fully neutral, 
e.g. in the dense parts of the interstellar medium, because 
ion-ion interactions, molecules and dust may then be impor- 
tant. 

We encourage the use of our website which pro- 
vides photo- ionised, non-equilibrium ionisation and cool- 
ing data as a function of temperature for a larger range 
of redshifts, meta l licitie s, and densities, including both the 
lHaardt &; Madaul (|200lh and (2012) backgrounds, than pre- 
sented here. We will continue to update this website with 
the release of new versions of CLOUD Y containing updated 
atomic data (e.g. iLvkins et all 12013). We also provide col- 
lisional and equilibrium results and provide some examples 
of how the calculations can be compared to observations. 

Diagnosing the physical state of the metal-enriched in- 
tergalactic and circumgalactic media where cooling times 
are short requires tracking non-equilibrium ionisation in the 
presence of the ubiquitous extra-galactic background. 
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APPENDIX A: EQUILIBRIUM IONISATION 
COMPARISONS WITH CLOUDY 

We compare equilibrium ionisation fractions calculated us- 
ing our method assuming a z = 1 HM01 background to 
those of CLOUDY ver. 10.00 for a range of diffuse gas den- 
sities and temperatures in Figure IA1I Red and black con- 
tours indicate ion fractions for our method and CLOUDY, 
respectively. Black contours overlap red contours in most 
cases, indicating excellent agreement. The shading indicates 
the ratio log [/rhisWork //cloudy] where either ion fraction 
is > 10 -3 if a metal ion or > 10 -8 if H or He; green indi- 
cates a good fit, while purple (orange) indicates our method 
yields too low (high) f Xi compared to CLOUDY. Note that 
our shading is very sensitive to differences between the two 
methods in order to highlight slight inconsistencies, and of- 
ten the contours nearly overlap when shading is off- green. 

For primordial species, the two methods agree well, ex- 
cept for Hi at T ^15,000 K and nu > 0.1 cm -3 where 
collisional ionisation causes a rapid transition in tempera- 
ture between neutral and highly ionised. We assume case A 
recombination in our method aimed to track diffuse ionised 
gases, and case B recombination is more applicable when the 
gas is optically thick to Lyman series photons. Nonetheless 
the two sets of contours still nearly overlap. 

Ovil and Nevin show excellent agreement, except at 
higher photo-ionised densities where f x . drops off rapidly 
with increasing density. This is also apparent in Figure [2] 
and is of little practical concern. C IV also shows excellent 
agreement, although slight deviations occur where collisions 
dominate below T = 10 5 K, which is also of little concern 
since /civ is very low here. Mgn shows good agreement, 
but reflects the disagreement with Hi at the same range of 
densities around T = 15,000 K, owing to spurious Hi and 
electron densities. Again, the contours show excellent agree- 
ment. The bottom three panels of Figure [ATI show a selec- 
tion of some of the worst agreeing species, which likely owe to 
more complex treatment of radiative and dielectronic recom- 
bination in CLOUDY ver. 10.00. These species are also less 
studied in laboratory experiments and the relevant atomic 
data is therefore often less well constrained. The excellent 
agreement we demonstrate between cooling efficiencies cal- 
culated by CLOUDY and us indicates that these ion abun- 
dance differences are not important for the dynamics of gas, 
but they may be important for the diagnostics of the gas 
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Figure Al. A selection of equilibrium abundances of common ions compared between our method (red contours) and CLOUDY (black 
contours). The colour bar shows the ratio of the ion fractions predicted by the two methods. In both cases, we turn on Auger ionisation 
and charge transfer. In general the agreement is very good as can be seen by comparing the black and red contours, although differences 
are highlighted by the very sensitive colour scaling that we choose. Where photo-ionisation is dominant for Hi and Hell, the agreement is 
excellent, but small differences occur in the regime where collisional ionisation dominates and Hi becomes primarily neutral; this has to 
do with our assumption of case A recombination, which breaks down in this regime. In general, Li-like and He-like ions (C IV, O VII, and 
Ne Vlll) show good agreement although slight differences can exist at higher photo-ionised densities where f x i drops off precipitously with 
increasing density. Mgll shows good agreement except where Hi fractions are underestimated, and the electron density is thus different 
between the two cases, which affects low ions like Mgll in the regime where collisional ionisation becomes important. We also show 
some species in the bottom row that are more spurious, likely owing to a more complex set of radiative and dielectronic recombination 
coefficients used in CLOUDY (Sin, Fen, and Fexvn). 
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if these particular species are observed. Comparisons with 
CLOUDY for additional ions can be found on the website 
: / / noneq.strw.leiden univ.nl| 
In Figure IA2I we demonstrate the effect of Auger ion- 
isation on O VI at z = 1. In the left panel we have Auger 
turned off both in CLOUDY and our method; the centre 
panel has Auger turned on in CLOUDY only; and the right 
panel shows it turned on in both. In the case of oxygen, 
photo-ionisation of the Is shell of O I through O V leads most 
often to the removal of two electrons, which is why O VI is 
more ionised with Auger ionisation turned on, especially at 
photo-ionised temperatures at densities above the peak of 
the photo-ionised O VI fraction. Hence, Auger ionisation is 
important for O VI strengths at these densities and temper- 
atures, and would be about 2x lower if Auger ionisation is 
ignored. 

We perform an analogous comparison for O I and the 
effect of charge transfer (CT) in Figure IA3I as we did for 
O VI and the Auger effect in Figure IA2I from left to right 
the panels show no CT in either case, CT only in CLOUDY, 
and CT in both. If CT is not turned on, the O I fraction is 
underestimated by an order of magnitude at temperatures 
where it is found. O I shows extreme sensitivity to charge 
transfer owing to O I having a very similar ionisation poten- 
tial as Hi. 
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Figure A2. Comparison of Ovi fractions predicted using our method calculated in equilibrium (red contours) and CLOUDY (black 
contours). Contours correspond to logarithmic ionisation fractions, and colour shading indicates the ratio of our method's ionisation 
fraction over that of CLOUDY, according to the logarithmic colour bar. From left to right, we show the Ovi fractions without Auger 
ionisation in either case (left), with Auger ionisation on for CLOUDY only (centre), and with Auger ionisation on in both methods 
(right). Auger ionisation boosts the Ovi fraction at photo-ionised temperatures at the high density end. 
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Figure A3. A similar plot as Figure lA2l but for Ol and the effect of charge exchange. Red contours indicate our ionisation fractions, 
and black contours indicate CLOUDY results. From left to right, we show the Ol fractions without charge transfer in either case (left), 
with charge transfer on for CLOUDY only (centre), and with charge transfer on in both cases (right). Charge transfer is very important 
for O I owing to its similar ionisation potential with H I, and it is also significant for many of the other lower ions. 
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